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! Abstract 
^ , 

O ! We consider the perturbative treatment of the minimally coupled, massless, self- 

^ \ interacting scalar field in Euclidean de Sitter space. Generalizing work of Rajara- 

^ ■ man, we obtain the dynamical mass m? oc y/XH^ of the scalar for non-vanishing 

PnJ ■ Lagrangian masses and the first perturbative quantum correction in the massless 

case. We develop the rules of a systematic perturbative expansion, which treats 
the zero-mode non-perturbatively, and goes in powers of V^. The infrared di- 
^ I vergences are self-regulated by the zero-mode dynamics. Thus, in Euclidean de 

Sitter space the interacting, massless scalar field is just as well-defined as the 
massive field. We then show that the dynamical mass can be recovered from 
the diagrammatic expansion of the self-energy and a consistent solution of the 
Schwinger-Dyson equation, but requires the summation of a divergent series of 
loop diagrams of arbitrarily high order. Finally, we note that the value of the 
long-wavelength mode two-point function in Euclidean de Sitter space agrees at 
leading order with the stochastic treatment in Lorentzian de Sitter space, in any 
number of dimensions. 



1 Introduction 



It is well-known that the free, massless, minimally coupled scalar field in the de Sitter 
background space-time cannot be defined in a de Sitter invariant way due to infrared 
(IR) divergences [Tj. This occurs in any number of dimensions ci, since the mode integral 
J d'^~^k that defines the Wightman two-point function always diverges logarithmically 
for |/c I — )■ 0. In the expanding cosmological coordinate frame the divergence arises from the 
red-shifting of modes, which leads to a pile-up of long-distance modes at late times. But 
non-interacting fields are not very interesting. The question arises whether upon turning 
on an interaction of the scalar quantum field, no matter how small, the IR problem could 
somehow cure itself. For this to happen, the interaction must effectively become non- 
perturbatively strong among the long-distance modes. If so, the non-perturbative dynamics 
may be too complicated to be solved with analytic methods. However, it may also be 
that after a suitable resummation or reorganization of the expansion in the interaction 
strength, the interacting, massless, minimally coupled scalar field lends itself to a well- 
defined, systematic treatment. 

Various previous results suggest that in scalar field theory with a quartic self-interaction 
— A0^/4! the originally massless field acquires a dynamical mass m'^y^ oc \fXH^, where H 
is the Hubble constant of de Sitter space-time, which indeed regularizes the IR divergence. 
Starobinsky and Yokoyama [2] treat the long-distance fluctuations of the field as a classical 
random field that satisfies a Langevin equation. The associated Fokker-Planck equation is 
solved for large times by a probability distribution that results in finite correlation func- 
tions. Another approach uses the Schwinger-Dyson equations and obtains the dynamical 
mass from a self-consistent solution. In the mean-field or the large- limit [3i,i4j the self- 
energy can be restricted to the one-loop, tadpole diagram. Garbrecht and Rigopoulos [5] 
analyzed the various in-in propagators in the CTP formalism and found that the large- iV 
result is modified by the two-loop self-energy, but remarkably, no further contribution arises 
beyond two loops due to systematic cancellations in the CTP index sums. However, while 
both formalisms agree on the parametric size \fXH^ of the dynamical mass squared, the 
two exact results from [2l[5] disagree on the numerical prefactor. Neither of the formalisms 
so far explains how to compute sub-leading terms systematically. 

The present work is motivated by the attempt to resolve the difference between the 
classical stochastic and diagrammatic/Schwinger-Dyson approach. For reasons that will 
become evident it is much simpler but still instructive to investigate the issue in Euclidean 
de Sitter space, which is simply the sphere S'^. In an elegant paper Rajaraman [6J consid- 
ered the functional integral on the sphere and identified the zero-mode integral as the origin 
of non-perturbative dynamics. He computed the two-point function of the zero mode, which 
can be related to the dynamical mass. In this paper we extend the functional-integral ap- 
proach and use 2PI methods to compute the exact self-energy, which is the central quantity 
in the diagrammatic approach. Our main results are as follows: 

• We formulate the rules for a well-defined perturbation expansion of correlation func- 
tions of the massless, minimally coupled scalar field in Euclidean de Sitter space. The 
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expansion parameter turns out to be a/A instead of the coupling A of the standard 
perturbation expansion. 



• For the massive scalar field we obtain the dependence of the dynamical mass on the 
Lagrangian mass m; for the massless field the leading 0{\/X) correction. 

• We show that the dynamical mass can be obtained from the loop expansion of the 
self-energy after summing a divergent series to all orders in the loop expansion. 

• It seems to have gone unnoticed that the Euclidean dynamical mass |:.6J agrees with the 
result from the stochastic approach. We show that this is true in an arbitrary num- 
ber of space-time dimensions despite the fact that the relevant dimension-dependent 
quantities are apparently unrelated. 

The interacting, massless, minimally coupled scalar field is therefore perfectly well-defined 
on the de Sitter background. For A ^ 1 there is a systematic weak-coupling expansion. 
The reason why this is possible despite the fact that the zero mode is truly strongly coupled 
is that the infrared theory consists of a single degree of freedom (the zero mode), whose 
dynamics can be solved exactly. What all this implies for Lorentzian de Sitter space is less 
clear. We must leave this important point to further investigation. 

2 Scalar field in Euclidean de Sitter space 

Euclidean de Sitter space is obtained from (i- dimensional de Sitter space in global coordi- 
nates with line element 

ds^ = dt^ - — cosh^(m) dnd-i (1) 
by defining i = —jj ~ and assuming 27ci/H periodicity in t, which turns ds"^ into 

ds^ = -g^^dx^dx" = -jp {dO^ + sin^ 9 dQa-i) ■ (2) 

Thus, Euclidean de Sitter space is equivalent to the (i- dimensional sphere with radius 1/H. 
Because the sphere is compact, functions admit a discrete mode expansion in spherical 
harmonics. In d dimensions the spherical harmonics (x) are labelled by the integer 
index vector L = {L, Lj^^i, . . . , Li) with L > L^-i > ■ ■ ■ > L2 > \Li\ and satisfy (see, 
e.g., 0) 

aY^ = ^ dAVgg'^'d, y^) = -h'l{l + - i)y^ (3) 

as well as the orthogonality relation 

J d'x./^)Yl{x)Y^,{x) = (4) 
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The volume of Euclidean de Sitter space is 



V, = I d'x^) = = (5) 

Here Yq denotes the lowest harmonic, which is constant. 

We consider the minimally coupled, real scalar field with Euclidean action 

S = J d'x^ I ^ 9^''d,<pd,<j> + ^m'<j>' + ^ 0' } 

L 

where the second line follows from the mode expansion 

L 

and the orthogonality relation (jl]). From the quadratic terms of ([6]) we deduce the free 
propagator 

^ ^ ^ Yr(x)Y;(x') 

with 



(4^)d/2 r^^) " ' V 2 ' 2 '2' 4 



\(d-\\ ra^ 

y = 2 {1 — cos6 cos6' — sinO sinO' w ■ w' ) and ~ / ( — 2 — ) ~ Jp ' 

Here ?/ is the invariant distance on the rf-sphere S"^, and w, w' are two unit vectors on the 
sub-sphere S'^^^ with solid angle element dQ^-i in The second line of (|8]) is indeed the 
de Sitter propagator in the Bunch-Davies vacuum [51IS] in imaginary time. 

The free propagator is ill-defined for m = 0. The leading term for small m is 

^ , ,x m-i-o H'^ r{d) 1 

Gfree X,X ^ —J- X — , 10 

(47r)"/"' r(|) m"' 

which, as can be seen from the first line of (|8]), originates only from the zero mode. Let us 
separate the constant zero mode from the field by defining 

0(X) =0O + 0(X). (11) 
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The free propagator is the sum of the zero mode and non-zero mode propagator, since 
cross terms vanish by angular momentum conservation. The free zero-mode propagator 
equals the right-hand side of ffTOl) for any m ^ , while the non- zero-mode propagator has 
a well-defined massless limit. In c? = 4, 

G,Ux,x') = ^ - llnf - 1^ (m = 0). (12) 



47r2 \y 2 4: 

From now on we consider the massless, scalar field, m = 0, unless mentioned otherwise, 
and assume A <^ 1. The free zero-mode propagator is not defined, which is not surprising, 
since the zero mode has no quadratic term in the action ([6]). The zero- mode must be 
treated non-perturbatively. Let 

G{x, x') = Go + G(x, x') = (0o0o) + (0(x)0(x )) (13) 

be the exact two-point function of the interacting theory. Since Go is constant, we may 
write 

Comparison with the L = term in the first line of ([8]) suggests that we identify ^dyn "with 
the dynamical mass of the originally massless scalar field, generated by the self-interaction. 
Note that this interpretation should be regarded with some caution, since the value of rri^y^ 
is not related to the decrease of correlation functions at large separation 1/mdyn- In fact, 
1/mdyn corresponds to distances parametrically larger than the radius of the sphere, which 
carry no meaning. Similarly, in Lorentzian de Sitter space a dynamical mass of order 
A^/^i/ is related to super-horizon correlations. Nevertheless, as will be seen below, a 
finite value of Go regularizes the IR divergence of the massless field and allows us to define 
a well-behaved perturbation expansion. 



3 Perturbation expansion on the sphere 

In [6] the zero-mode two-point function (0o0o) was computed by evaluating the dominant 
contribution to its functional-integral representation, which gives 



m 
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V4! r(|) V 



VXH''/^ (15) 



In the following we generalize this approach. We show that both, (0o0o) and (0(a;)0(x')) 
have well-defined perturbation expansions in y/X, and provide a set of Feynman rules for 
this expansion. 

The generating functional is conveniently written in terms of two separate sources Jq, 
J{x), for the zero- and non-zero-mode field, respectively: 

Z[Jo, J] = N j V[(Po]V[^] exp (^-S - j d^x^i^M^ + ^0)) 
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exp ( —Sj 



'int 



Zo[Jo]^frec[J], (16) 



where N is defined such that 2'[0, 0] = 1. Here 



^int[</'o,0] = ^ j d'x^(^^^ + 4<j,o<j^' + 6<f)l<f)''j , (17) 

and the term proportional to (p^cj) vanishes since / d'^x^Y^ix) = for L > 0. The key 
point is that the term ^ 0q is not included in Sint, but must be part of ^ol-^o], since in 
the absence of a mass term for the scalar field the quadratic term in the zero-mode action 
vanishes, and the integral over (pQ in Zo[<^o] does not converge for large field values [B]. 
Hence, 

Zo[Jo] = Noj V[<Po] exp j d'^x^ (^^ 0^ + Jo0o)) , (18) 
while the generating functional for the free non-zero-mode field, 

ZireeiJ] = N j V[$\ CXp j d^X^ (^9^'' O^^d,^ + J 

= expQ j d'^xy/g j d'^y^ J{x)Giree{x,y)J{y)^ , (19) 

is a standard Gaussian functional integral. The zero-mode functional integral is simply an 
ordinary one-dimensional integral. Moreover, 0o and Jo are independent of x, so j d'^x^/g = 
Vd in f|T8|) . The integral over 0o can be evaluated exactly. Introducing 

A = ^Vd, Jo = JoVd, (20) 

we find 

Zo[Jo] = No I 6/00 exp f - Jo</>o 

oo 

2' 4' 256A ; 2r Q) ^ 1^4' 2' 256A J ' ^ ' 

where 0-^2 denotes a hypergeometric function. One easily checks that 

H VA r (4) 

reproduces (|T5l) as is should be. Similarly, (0o")o follows from taking the appropriate 
number of derivatives. The index on the bracket means that the computation is done 
with the zero-mode functional Zq alone. The full zero-mode n-point functions computed 
from Z in (ITB]) receive sub-leading corrections, as discussed below. 
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It is now straightforward to develop a systematic perturbative expansion of fll6l) and 
the corresponding Feynman rules. From (ITS]) it follows that every zero-mode fields counts 
as A~^/^, while has the standard counting 1. The interaction ^^o*/*^ therefore counts as 
\/A. In general, since there is always an even number of 0o involved, correlation functions 
have an expansion in a/A. The rules are as follows: For a given correlation function expand 
f lT6|) to the desired order in A using the above counting rules. Perform the standard Wick 
contractions of pairs of non-zero mode fields. This can be represented in terms of lines 
and vertices in the usual way. However, no Wick contractions are to be performed for the 
zero-mode fields. Instead, collect all factors of (pQ and compute the expectation value of 
exactly. 

As an example, we evaluate the first correction to the zero-mode and non-zero-mode 
two-point functions. For the zero-mode case, we have 



(0O0o) = (0o)o + mo X 

= (0o)o - ^ (0o)o I d'x^GlZi^,x) + 0{^/X), (23) 

where the black square represents the ^(pl(f)'^ vertex. If we define {(pocpo) = l/(Vdm^yj^) as 
before in f lT^ and denote the previously obtained leading-order expression f|T5|) by m^y^ q, 
the previous equation translates into 



^dyn ~ '^dyn.O 



1 + v^-il^il- ^,G}Z{x,x) + 0{\) \ , (24) 



where (0o)o = 1/(4A) has been used. While the leading expression m\y^Q is unambiguous, 
the first correction depends on the UV subtraction that defines the coincident non-zero- 
mode propagator G\^^^{x^x). To be specific, consider the case of four dimensions. In 
dimensional regularization one needs to take x' — )■ x before expanding the d = 4 — 2e 
dimensional propagator Giree{x-,x') around c? = 4, in which case 
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? + 21n-^-7 + 47s 



(25) 



The MS-renormalized value corresponds to this expression with the pole term in 4 = 
7 — 7b + ln(47r) in brackets subtracted. Further, /x is the MS renormalization scale. We 
note that the dynamical mass is not by itself a physical quantity. While the leading term is 
unambiguous, the first correction involving the propagation of non-zero modes is scheme- 
and scale-dependent. 

The non- zero-mode two-point function is the free propagator in leading order. Including 
the C(V^) correction it reads 

(0(x)0(x')) = Gfree(a;,a;') + (0o)o X ■ 
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Giree{x,x') - ^ (0o)o / Z ^/g G iree{,X , z)Giree{,Z , x') + 0(A) 



H^L{L + d-r 



1 - 



2H^L{L + d-r 



+ 0{X) 



(26) 



The expression for the leading correction differs from The diagram computed there is 
part of the sub-leading 0{X) correction. In four dimensions the leading term equals ( IT^ . 
and the leading \/X correction can also be summed to give 



A(0g)o 
487r2 



!.,(.)+ (--2jln(l-.)- - + - 



(27) 



where u = 1 — y/A. 



4 Schwinger-Dyson equation 

We now return to the approaches pursued in [3-5J which are based on evaluations of the 
scalar-field self-energy and the Schwinger-Dyson equation 

D^G{x, x') = i (5^^^) (x - x') + [ d'^w>/g(w)U{x,w)G{w,x'). (28) 

We project this equation on the zero-mode component by integrating over x' and using the 
identity 

d'xyW)Yd^) = ^^, (29) 
which follows from (HI). This results in 



= 1 + v;,'noGo. 



(30) 



In the spirit of the 2PI formalism (see below) we regard the self-energy as a functional of 
the exact propagator and derive it from the functional derivative 



U(x,x') 



rprest 
'^^ 2PI 

6G{x, x') 



(31) 



of the 2PI effective action with the classical and one-loop term subtracted (as denoted by 
"rest") lim. The loop expansion of the effective action is given by 



prest 
2PI 






J_ 

48 



48 




C(A^ 



1 

128 



32 



oo „ „ 

Y/x^{-XT j d'^x^^/^... j rfW^ G(.,.). . .G(.,.)^ , 



2n factors 



(32) 



where a„ denotes the combinatorial factor associated with a diagram, as given in the first 
hne, and where we have given the exphcit diagrammatic representation up to the five-loop 
order. Instead of appealing to the 2PI formalism we could have written down the self- 
energy diagrams directly with the proviso that all internal lines are exact rather than free 
propagators. 

The power-counting rules of the previous section tell us that the leading contribution 
to Fgpj and IIo is obtained from pure zero-mode diagrams, that is, every full propagator is 
replaced by Gq. Since every loop brings one factor of A from the new vertex and adds two 
propagators, which each count as 1/ a/A, we conclude that every order in the loop expansion 
contributes to the leading term. This shows that the loop expansion must be summed to all 
orders to obtain the correct value of the "dynamical mass" in Euclidean de Sitter space. 

To see this explicitly, note that upon plugging the derivative of ( 15^ into ( 15U]) . the latter 
equation can be solved for Gq, which yields the value of rri^y^. (More precisely, m'^y^ Q since 
by replacing G by Go in ^ 132]) we pick up the leading term only.) The functional derivative in 
f l3T|) eliminates two integrations in f l32|) such that an ra-loop self-energy diagram contributes 
an(-A)'^V7^^Gg""^ to dSnDEl Therefore, with the ansatz Go = (XVd)-^/^ x z, where z is 
a number to be determined, and the explicit expression (15^ for the loop expansion up to 
the four-loop self-energy, the Schwinger- Dyson equation fl5U]) turns into 

z^ z"" 

= l-y + ^-^ + ^ + . ..-/(.), (33) 

where the term 2;^" corresponds to the sum of n-loop self-energy diagrams. In terms of z, 
the "dynamical mass" is given by 



y/X d=4 1 VOXH^ 
~ ~z^ 47r 



Keeping only the one-loop tadpole diagram in (133|) . which corresponds to truncation after 
the quadratic term, we obtain z = a/2 and m"^^^ = \/3XH'^/ (Att) (in d = 4), which coincides 
with the mean- field result [2] and the one- loop result in Lorentzian de Sitter space [3H5]. 
In the two-loop approximation, the quadratic equation for z"^ does not yield a real solution 
for z. Thus, at two loops there is a difference between the solution of the Schwinger-Dyson 
equation in Euclidean de Sitter space and the solution to the corresponding equations for 
the closed-time-path propagators in Lorentzian de Sitter space 0. Continuing to higher 
orders in f l33|l . at three loops, we find z = 1.166 . . ., while at four loops (which includes the 
last term shown explicitly in f l5^ ) there is again no solution. Thus, the loop expansion 
does not seem to converge to the exact value ffT^ . which corresponds to 

r 

4!— )4v = 1-65580.... (35) 



-'^Alternatively, we can substitute G Go in the expression ((32|) for Fj'pj and define the zero-mode 
self-energy as the ordinary derivative 2dT^^l/dGo. With this convention Hq as defined above must be 
replaced V^Uq and the Schwinger-Dyson equation (j5D|) takes the form = 1 + IIoGo. This convention is 
adopted in Section 6 below. 
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The question arises how the exact result that is obtained easily from the functional integral 
is recovered diagrammatically. Clearly, we need the expansion of f{z) to all orders. But 
this cannot be obtained from fl32|) . While the integrations are trivial in the zero-mode 
approximation, the diagram topologies and computation of combinatorial factors become 
too comphcated. 



5 Zero-mode dynamics in the 2PI formalism 

In the following we exploit the 2PI formalism [Tm[TT| to derive an expression that generates 
the perturbative expansion of the zero-mode self-energy to any desired order. We focus 
on the zero-mode dynamics which alone is responsible for the leading contributions as 
discussed above, and hence set to zero. In this section we drop the subscript "0" , since 
all quantities are understood to refer to the zero mode. 
The generating "functional" in the 2PI formalism is 



Z[J,R] = e-^[^'^] = N f exp (-S[4>] - j J{x)^{x) - - f 



R{x,y)<^{x)<^{y) 



x,y 

oo / 



— oo 



/ dcj) exp ( -A0^ - J0 - -Rcj)"^ ) , (36) 



where A and J (equal to Jo) have been defined in (120|) . [...] = J d'^x^J g{x) [...], and 

R = V^R. (37) 

The simple one-dimensional integral in the second line of (136|) applies since the zero-mode 
field is constant. The exact propagator in the presence of the external sources can be found 
from the relation 

^ = \ {^l + G) , (38) 

where 0ci = dW/dJ is the field expectation value. Since the functional integral is an ordi- 
nary one- dimensional integral, the functional derivatives are actually ordinary derivatives. 
It follows from (pT]) that the field expectation value vanishes in the absence of the source 
J, that is, the symmetry —0 is not spontaneously broken. This remains true for 
J = and R ^ 0. Since eventually we are interested in the theory in the absence of 
external sources we now put J = and consequently 0ci = 0. In this case we find the 
closed expression 




Z[0, R] = e--[°^«l = ^^^^^^ exp I I ir, I f;^ I , (39) 
where K^{x) denotes the modified Bessel function of the second kind. From (138|) we obtain 
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0.1 0.5 1.0 5.0 10.0 50.0 100.0 



Figure 1: Dynamical mass of the zero mode as function of the Lagrangian mass. 

The previous equation gives the exact zero-mode propagator or, equivalently, the "dy- 
namical mass" of the zero mode in the presence of the propagator source R. We note 
that since R is constant, V^R has an equivalent interpretation as a Lagrangian mass 
for the scalar field. Hence (l4Up provides the "dynamical mass" of the scalar field for ar- 
bitrary m, generalizing the expression f|T5|l for the massless case [6]. The dependence of 
m^yjj on is sketched in Figure [U At large the "dynamical mass" asymptotes to 

with ordinary perturbative corrections of order A, as should be expected, since the 
infrared enhancement that renders the zero-mode dynamics non-perturbative is cut off for 
a sufficiently massive scalar field. For — )■ 0, the "dynamical mass" tends to the value 
f[T^ . The pre-asymptotic corrections can be determined easily by expanding f HU]) around 
the corresponding limits. 

6 Diagrammatic zero-mode self-energy to all orders 

We now determine the self-energy that is needed to solve the Schwinger- Dyson equation. 
In the 2PI formalism, the Schwinger- Dyson equation reads 

G-\R) = Gi^l + R-U{G). (41) 

The inverse free propagator is G^^^ = — 12A0^i, which vanishes for J = 0. This refiects 
once more the fact that the free propagator of the massless scalar field is ill-defined in the 
absence of external sources. Thus, 

U{G) = -^ + R{G). (42) 

It follows that the diagrammatic expansion of the self-energy is obtained by inverting G{R) 
given in (HUI) . and expanding it in powers of A. 
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n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


an 


1 

2 


1 

6 


1 

4 


5 
8 


101 
48 


279 
32 


8143 
192 


271217 
1152 


374755 
256 


51151939 
512 



Table 1: The first ten series coefficients a^. 

While a closed expression for the inverse of G{R) may not exist, we can solve for the 
expansion in A by making the ansatz 

R{G) = ^ + a^WdG + a^X^V^G^ + . . . , (43) 
G 

where the ffist term is required by f H2]) to obtain a regular perturbative expansion of 
n(G). The term a„A"'VJ*G^"'"'"^ represents the sum of the n-loop diagrams to the zero- 
mode self-energy, expressed in terms of the exact zero-mode propagator. The definition of 
the expansion coefficients a„ is chosen such that with the definition (134|) of z the function 
f{z) defined in (133|) is given by 

oo 

/(^) = l-^a„^2n_ (44) 

n=l 

Plugging the ansatz f H5]) into f HU]) and matching coefficients in the expansion in A (equiv- 
alently, in G), we find the a„. The ffist ten terms are shown in Table [TJ We note that 
the first four agree with ( 133|) obtained from the combinatorial factors of the lowest-order 
Feynman- diagram topologies. We determined the exact coefficients up to n = 200, which 
turn out to be rational numbers of increasing length. 

Inspection of the a„ coefficients shows that they form a sign-alternating, factorially 
divergent series with 



n—¥oo 



0.201... X ( -^"j + ...] . (45) 



The divergent behaviour arises because the expansion in small G corresponds to an expan- 
sion of G{R) around R = oo, see fH3|l . while the value of m^y^ is related to G{R) at R = 0. 
The divergent series is also the reason why we did not obtain a reasonable approximation 
to wijyj^ from the low-order approximations to the self-energy. 

It remains to show that the Schwinger- Dyson approach is consistent with the exact 
result f[T^ for the "dynamical mass" , which requires summing the divergent series. To this 
end we construct the Borel transform of f{z) 



!?[/](«) = 1 a„-, (46) 

n=l 

such that the Borel sum F{z) of f{z) is given by 

1 

F{z) = - due-^l^' B[f]{u). (47) 
^ Jo 
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n 


3 


6 


9 


12 


24 


48 


96 


z^: from F 


1.65709 


1.65635 


1.65581 


1.65580 


1.65580 


1.65580 


1.65580 


from / 


1.71012 


1.66262 


1.65723 


1.65618 


1.65581 


1.65580 


1.65580 



Table 2: Solution of F{z^) = (second line) and /(z*) = (third line) employing a diagonal 
{n,n) Pade approximation constructed from the first 2n series coefficients a^. The exact 
value is 2, = V^r (|) /T (|) = 1.65580 .... 

Since the series is sign-alternating, we expect B[f]{u) to exhibit a singularity on the nega- 
tive axis, but without a closed expression we do not know the precise singularity structure 
of the Borel transform of /. Given f l40p . it is reasonable to assume that it is analytic in a 
vicinity of the positive real axis such that the Borel integral is well-defined, and to assume 
that the Borel sum F{z) equals the original function f{z). 

Term- by-term integration of the series expansion of B[f] (u) simply returns the divergent 
series expansion of /. We therefore resort to a standard trick [12] and construct a Pade 
approximation from the truncated series expansion. More precisely, we use the first 2n 
coefficients of the expansion of B[f], not counting the "1" in (jM]), and construct the 
diagonal (n, n) Pade approximant. We use this approximation to B[f]{u) in the Borel 
integral ( H7|) and obtain F{z) by numerical integration. We then solve the equation F{z^,) = 
0, see (133|) . to determine the "dynamical mass". Alternatively, we can determine the 
solution of f{z^) = from Pade approximants to the expansion (Hlj) of f{z) directly, 
without going through the Borel transform. The results are shown in Tabled The solutions 
are seen to quickly approach the exact result ( 135|1 . especially when the Pade approximation 
is applied to the Borel transform, in which case one-permille accuracy is reached already 
for n = 3. This demonstrates that the diagrammatic approach via the 2PI Schwinger-Dyson 
equation reproduces the path-integral result, as it must be, but only after summation of a 
divergent series expansion to all orders. 

7 Stochastic approach in d dimensions 

The methods applied above do not extend to Lorentzian de Sitter space, which is non- 
compact, and does not allow to identify the (leading) infrared dynamics with the one of 
a single zero-mode degree of freedom. However, quite some time ago Starobinsky and 
Yokoyama [2j suggested that the long-wavelength part of the scalar field can be treated 
as a classical stochastic variable, which satisfies a Langevin equation with a random force 
provided by the short-wavelength modes. Here we show that this leads to the same value 
for the two-point function of the long- wavelength field as the zero-mode two-point function 
in Euclidean de Sitter space, in any number of dimensions d. This intriguing coincidence 
seems not to have been noted before. 

Following [2| we divide the scalar field into (j) = (po + (f), where (po contains all long- 
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wavelength modes with wave number k < eaH with a the scale factor and e ^ 1 the 
parameter that separates long from short wave-lengths. From the field equation it follows 
that 00 satisfies the Langevin equation 

Mt,x) = -^^-L_ V'(0o) + f{t,x). (48) 

Here V{(j)) is the scalar field potential and / the stochastic force 

f = } = eaH^ J -^-^ 6{k - eaH) (akMt)e''-'' + cc) (49) 

generated by the short- distance modes. At leading order, we can neglect the self-interaction 
of the short- distance modes. The fluctuations satisfy {f{ti,x)f{t2,x)) = a6{ti — ^2) with 



(27r)'^-ir(^^) 71 ^ ' 

The first factor arises from the volume of the d — 2 dimensional momentum shell k = eaH, 
the second from the long-wavelength limit (since k = eaH ^ aH) of the Bunch-Davies 
mode functions (f)k{t)- The Fokker-Planck equation for the one-particle probability density 
P[(p] associated with ( HHj) is 

dP Id a d'^ 

which admits the stationary late-time solution 

in terms of which the two-point function of the (constant) long-wavelength field is given 
by 

{^M = N £ ^= exp {-jj^ V(^)) . (53) 

This precisely agrees with fl22l) (for V{(f)) = ^(p^) provided the dissipation and fluctua- 
tion coefficients in the Fokker-Planck equation are related to the volume of li-dimensional 
Euclidean de Sitter space with radius 1/H hy 

id-l)Hx^ = ^, (54) 

which can be easily verified. Hence, the long-wavelength two-point functions (and therefore 
"dynamical masses") are the same, as claimed. One may wonder why the result could be 
derived from zero-mode dynamics alone in Euclidean de Sitter space, while the fluctuations 
originated from the short-wavelength modes in the stochastic approach. However, the result 
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from the latter is independent of e in the above approximation, as it should be, and the 
stochastic force is generated by wave numbers k = eaH. We can take e arbitrarily small 
and conclude that in the leading approximation the main contribution to the stochastic 
force can be assumed to originate from the boundary between long- and short wavelengths, 
which can be taken to be deep in the infrared. 

We note that the stochastic approach can be derived rigorously from the full quantum 
dynamics in the leading logarithmic approximation in a [13] , which is equivalent to keeping 
the leading infrared-enhanced terms in Euclidean de Sitter space. But unlike the Euclidean 
case discussed in the present paper, a systematic method for calculating corrections around 
the Lorentzian result of [2j is not known. 

8 Conclusion 

In this paper we considered the perturbative treatment of the minimally coupled, massless, 
self-interacting scalar field in Euclidean de Sitter space. Generalizing the work of Rajara- 
man [6], we obtained the dynamical mass oc \/XH^ of the scalar for non- vanishing 
Lagrangian masses, as well as the first perturbative quantum correction in the massless 
case, and developed the rules of a systematic perturbative expansion, which after treat- 
ing the zero-mode non-perturbatively, goes in powers of V^. We then showed how the 
dynamical mass can be recovered from the summation of the diagrammatic expansion of 
the self-energy and a consistent solution of the Schwinger-Dyson equation. This clarifies 
the relation between the path-integral and diagrammatic treatment, and implies that so- 
lutions based on truncations of the loop expansion can at best be approximate. With the 
proper exact treatment of the zero mode, the reorganized perturbative expansion is free 
from infrared divergences, which are present for the free, minimally coupled scalar field 
in Euclidean de Sitter space. The interacting, massless field is therefore well-defined, and 
the rules for generating the systematic perturbative expansion are almost as simple as the 
standard rules for the massive case. 

What this implies for Lorentzian de Sitter space is much less clear. We showed that 
the long-wavelength mode two-point function computed in the stochastic approach of [2] 
coincides with the exact Euclidean result in leading order in the expansion in v^. This 
strongly suggests to us that the dynamical mass of the self-interacting scalar field in de 
Sitter space can be obtained by some sort of analytic continuation from the Euclidean, up to 
higher-order corrections. It would be very interesting to derive this result diagrammatically, 
in the spirit of |5], and to understand how to develop a systematic expansion in Lorentzian 
de Sitter space. 
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